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ABSTRACT 

Fnll-polarization MERLIN and VLBA observations of the 1720-MHz maser 
emission from the OH molecnle in the snpernova remnant W44 are presented. 
Zeeman splitting has been directly resolved between right- and left-circnlarly 
polarized spectra indicating magnetic held strengths of |B| 0.50 mG. The 

position angle of the linear polarization is in good agreement with the orientation 
of the shocked gas at the snpernova remnant/molecnlar clond interface. Based on 
MERLIN data having a resolntion of 200 milliarcseconds (mas) and VLBA data 
with a resolntion of 15 mas, the masers are measnred to have deconvolved angnlar 
sizes of 50 to 350 mas (150 to 1000 AU) with compact cores 20 mas (60 AU) in 
size, consistent with theoretical expectation and previons observations. 

Subject headings: ISM: individnal (W44)—masers—snpernova remnants 


1. Introduction 


Spectral-line emission from the OH (1720 MHz) satellite transition (^n 3 / 2 , ^ = |,F = 
2—i>l) in snpernova remnants (SNR’s) was first observed by Goss (1968) toward W28 and 
W44. The SNR OH emission in the 1720 MHz line is non-thermal and is accompanied by 
absorption in the 1612-, 1665-, and 1667-MHz ground state hyperhne transitions of OH (Goss 
& Robinson 1968; Ball & Staelin 1968; Turner 1969; Robinson, Goss, & Manchester 1970; 
Goss, Gaswell, & Robinson 1971). Beginning with observations of W28 by Frail, Goss, and 
Slysh (1994), subsequent studies and surveys (Yusef-Zadeh et ah 2000; Yusef-Zadeh et ah 
1999b; Koralesky et ah 1998; Green et ah 1997; Frail et ah 1996; Yusef-Zadeh et ah 1996) 
have observed ~ 200 Galactic supernova remnants and found 22 remnants with associated 
OH (1720 MHz) maser emission. 

Elitzur (1976) suggests that the level inversion of these 1720 MHz SNR masers is due 
to collisional excitation. The physical conditions suggested for the collisional pump (n ~ 
10®cm“^, T ^ 90 K, Yqh ~ 10^®“^^cm“^) are observed to occur where a G-type (non- 
dissociative) shocked interface between the SNR and an adjacent molecular cloud is moving 
transversely across the sky {e.g. Elitzur 1976; Wardle 1999; Lockett, Gauthier, & Elitzur 
1999; see also Frail et al. 1996). The temperature, density, and OH column of the OH 
(1720 MHz) collisional pump are well established from observation [e.g. Glaussen et al. 
1997, hereafter G97; Hoffman et al. 2003a, hereafter H03). Furthermore, both Glaussen et 
al. (2002) and H03 hnd the intrinsic sizes of the masers to be consistent with predictions 
based on collisional pump theory (~10^^ cm, Lockett et al. 1999). Also, observations of the 
shocked emission from GO (J = 3 —2) near the OH masers in W28 and W44 by Frail & 
Mitchell (1998) and Arikawa et al. (1999) indicate that the masers originate in the post-shock 
gas, consistent with the collisional-pump model (see also Reynoso & Mangum 2000; Dubner 
et al. 2004; Lazendic et al. 2004). 

OH (1720 MHz) SNR masers exhibit relatively large angular size (approximately 100 mas, 
H03), modest flux density (< 1 Jy, e.g. Frail et al. 1996), and close relative angular spacing 
(> 100 mas, H03). The masers also exhibit ~ 10% linear and circular polarization [e.g. 
G97, H03). However, several masers may be “blended” into a single image feature in ob¬ 
servations with angular resolution > 500 mas, prohibiting the determination of the intrinsic 
polarization state of the maser emission. Observations using MERLIN (approximately 100- 
milliarcsecond [mas] beam) and the VLBA (approximately 10 mas beam) are required in 
order to determine the intrinsic emission properties of OH (1720 MHz) SNR masers. 

In a recent paper (Hoffman et al. 2005, hereafter Paper I), we presented full-polarization 
MERLIN and VLBA observations of approximately 30 OH (1720 MHz) SNR masers in W28. 
The observations of W28 in Paper I were shown to be relatively free from confusion due to 



spatial or spectral blending. Thus, the intrinsic polarization properties of the maser emis¬ 
sion were determined, indicating that: (1) the magnetic held orientation in the plane of the 
sky indicated by the linear polarization state of the maser emission is in good agreement 
with the magnetic held traced by synchrotron emission the SNR limb and (2) direct obser¬ 
vation of the Zeeman splitting between the right- and left-circularly polarized maser lines 
in OH (1720 MHz) SNR masers is a reliable indicator of the strength of the magnetic held. 
However, it should be noted that although observations using MERLIN and the VLBA pro¬ 
vide a comprehensive empirical understanding of the maser emission, the interpretation of 
the results using existing theoretical models is relatively unconstrained {e.g. Paper I; Gray 
2003). In particular, the theoretical derivation of the angle 6 between the line of sight and 
the magnetic held in the maser gas is not well constrained. 

The W44 SNR also contains a large number of masers (>25) compared with most OH- 
maser SNR’s, which typically contain only a few masers {e.g. Brogan et al. 2000, hereafter 
BFGT); an observational study of the W44 masers should yield a statistically signihcant 
sampling of OH (1720 MHz) SNR maser polarization properties for comparison with the 
conclusions of Paper I. W44 is a SNR of diameter approximately 26 pc in the radiative 
phase of evolution {e.g. Ghevalier 1999). The remnant has an age of about 2 x 10"^ yr and 
a distance of 3.1 kpc determined from observations of the interstellar absorption against 
the SNR {e.g. Radhakrishnan et al. 1972; Gaswell et al. 1975) and of the dispersion due 
to the ionized interstellar medium in the direction of the associated pulsar PSR B1853-I-01 
(Wolszczan et al. 1991). The molecular cloud with which the SNR is interacting has been 
observed in both H I and molecular lines {e.g. Wootten et al. 1977; Sato 1986; Koo & Heiles 
1995; Seta et al. 2004). The molecular cloud/SNR interaction in W44 appears to have a 
relatively simple geometry; the molecular cloud lies almost directly between the observer 
and the SNR (Kodaira et al. 1977; Wootten 1977). Seta et al. (2004) suggest that the 
SNR/molecular cloud interface has corrugations ~ 1 pc in length so a transverse shock 
front (suggested to be necessary for the observation of the OH masers) is present at several 
locations interior to the SNR limb (see their Fig. 11). Indeed, Figure 1 shows an image of 
the synchrotron radiation from W44 with the positions of the OH (1720 MHz) maser regions 
overlaid. The masers are observed across most of the face of W44 in contrast to the case 
of W28 (G97) where the masers occur along the limb bordering the molecular cloud at the 
northeast of the SNR. Furthermore, the velocities of the masers in W44 he within the range 
of approximately 3 km s“^ (a few linewidths) in contrast to the approximately 11 km s“^ 
(> 10 linewidths) over which the W28 masers are observed. The geometry and velocity 
structure of the W44 SNR/molecular cloud interaction appears to have been shaped by the 
creation of a wind-blown cavity prior to the supernova explosion {e.g. Koo & Heiles 1995). 

The OH (1720 MHz) SNR masers in W44 were previously observed by G97 at 1" resolu- 
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tion in an extensive (25 masers), full-polarization experiment using the VLA and by Claussen 
et al. (1999, hereafter C99). The maser emission was also observed in a comparatively small- 
scale (2 masers), circular-polarization experiment using MERLIN and the VLBA. In this 
paper, we continue the investigation of the polarization of OH (1720 MHz) SNR masers with 
full-polarization observations of most of the OH masers in W44 using the Multi-element Ra¬ 
dio Linked Interferometry Network^ (MERLIN) and the Very Long Baseline Array (VLBA) 
of the NRAO^. This paper will use the source numbering convention of Claussen et ah (C97 
Table 1). 


2. Observations 
2.1. MERLIN+Lovell 

The W44 masers were observed on 26 and 27 January and 28 February 2002 using 
the MERLIN radio telescope of Jodrell Bank Observatory for a total of ~ 10 hours. The 
observations used two array tracking positions: (1) the E and F maser region and (2) the B 
and C maser region (C97); neither the A or D region masers were observed using MERLIN 
due to time constraints. The observations were centered on a Doppler velocity of ulsr = 
45.5 km s“^ for a line rest frequency of 1720.52998 MHz. Seven antennas were used; the 
Mark H and Lovell telescopes at Jodrell Bank, the 32 m antenna at Cambridge, and the 25 m 
dishes at Knockin, Darnhall, Tab ley, and Defford. The very short (< 500 m) Mark H-Lovell 
baseline was not used. The baseline lengths of MERLIN range from 11 km to 217 km; the 
array is not sensitive to angular scales larger than 373; therefore no continuum radiation from 
the SNR was detected. The antennas have right- (R) and left- (L) circularly polarized feeds 
from which RR, LL, RL, & LR cross-correlations were formed. The correlator produced 256 
spectral channels across 250 kHz (44 km s“^). The spectra were Hanning-smoothed off-line 
yielding a velocity resolution of 0.36 km s“^. The visibilities were integrated for 8.0 s. 

The absolute amplitude calibration is based on observations of 3C 286. The bandpasses 
were calibrated based on observations of 3C 84. The phases were calibrated using frequent 
observations (with a bandwidth of 14 MHz) of the source 1904-1-013 (J1907-I-0127), which has 
a flux density of 0.3 Jy and a positional uncertainty of 8 mas. The uncertainty in the absolute 


^MERLIN is operated as a National Facility by the University of Manchester, Jodrell Bank Observatory, 
on behalf of Particle Physics and Astronomy Research Council (PPARC). 

^The National Radio Astronomy Observatory (NRAO) is a facility of the National Science Foundation 
operated under a cooperative agreement by Associated Universities, Inc. 
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positions of the masers is about 10 mas (depending on the signal-to-noise of the individual 
detections). The visibility phases and amplitudes were then (self-) calibrated using images of 
the Ell maser (approximately 3 Jy, C97). 3C 84 was observed at both 14 MHz and 256 kHz 
bandwidth allowing bandwidth-dependent calibration to be transferred. The position angle 
of the linear polarization response of the antennas was determined from observations of 
3C 286 and 1904-1-013. Linear-polarization calibration sources were not observed for the B 
and C region masers. The RMS noise in a line-free, single channel image is 20 mJy beam“^, 
in agreement with expected instrumental behavior. The FWHM synthesized beam of the 
images is 240 x 141 mas at a position angle 24°. 


2.2. VLBA+Yl 

The W44 OH (1720 MHz) masers were observed on 13 and 18 April and 9 and 15 
May 2003 with the ten antennas of the VLBA plus one antenna of the VLA for a total of 
approximately 18 hours. The antennas have right- (R) and left- (L) circularly polarized feeds 
from which RR, LL, RL, & LR cross-correlations were formed. The correlator produced 128 
spectral channels across a 62.5 kHz (11 km s“^) band yielding a 0.09 km s“^ velocity channel 
spacing. The observing band was centered on ulsr = 45.425 km s“^. The spectra were 
off-line Hanning-smoothed yielding a velocity resolution of 0.18 km s“^. The visibilities were 
integrated for 7.3 s. 

The amplitude scale is set using online system temperature monitoring and a priori 
antenna gain measurements. Bandpass responses and station delays were found from obser¬ 
vations of 3C 345. The phases were calibrated using frequent observations (fast switching) of 
the source J1838-I-0404, which is approximately 5° in angular distance from the masers with 
a flux density of 0.3 Jy and with a positional uncertainty of 3 mas. The uncertainty in the 
absolute positions of the masers is approximately 3 mas (depending on the signal-to-noise of 
the individual detections). The visibility phases and amplitudes were then (self-) calibrated 
using images of the Ell maser (6.6-Jy flux density observed with the VLA by C97; 1.2-Jy 
flux density observed with the VLBA by C99, see also §§3.1, 4.1). The position angle of the 
linear polarization response of the antennas is determined from observations of 3C 286 and 
J1751+0939. 

The baseline lengths of the VLBAj-Yl array range from 52 km to 8611 km; the array 
is not sensitive to angular scales larger than 0770. The correlated held of view with respect 
to the phase-tracking position is 15". Thus, in order to observe the W44 masers, two array¬ 
pointing positions and hve correlator-tracking positions were used, summarized in Table 1. 
The A region maser was not observed with the VLBA due to time constraints. Also, the data 
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were significantly contaminated with radio-frequency interference, especially at the Brewster 
and Saint Croix stations, causing decreased angular resolution in the northwest-southeast 
direction and decreased image sensitivity. The severe anisotropy of the angular response of 
the array was partially remedied through weighting (tapering) of the visibility amplitudes in 
the M,n-plane {e.g. Briggs, Schwab, & Sramek 1999). The synthesized beam of the resulting 
images is 27 x 13 mas at a position angle of —23°. The RMS noise in the hnal single-channel 
images is 25 mJy beam“^, in agreement with expected instrumental behavior. 


3. Results 

3.1. Angular Structure 

All calibration and imaging was performed with the AIPS^ software package. Figure 2 
shows the MERLIN image of the Ell maser region. Figure 3 shows the VLBA image of the 
F24 maser region. Tables 2 and 3 list the deconvolved angular sizes and total flux densities 
observed using MERLIN and the VLBA. 

All of the five masers observed using MERLIN for which the peak intensity in VLA 
observations (C97) is greater than 1 Jy beam“^ were detected (B7, Ell, E20, F23, and F24). 
In addition, the C8 maser (700 mJy beam“^ with the VLA) was detected in the MERLIN 
data. This flux density limited detection rate is consistent with the MERLIN observations 
of W28 (Paper I). The C8, Ell, and E20 masers have multiple angular components in the 
MERLIN images {e.g. Fig. 2), consistent with the observations of the Ell maser by C99 
using MERLIN. Similarly in W28 (Paper I), six masers in the MERLIN data were shown to 
have multiple components with relative angular separations < 1". The E20 and F23 masers 
were not detected in the MERLIN observations by C99, as discussed in §4.1. 

The current observations using the VLBA yield detections of two of the ten masers 
detected using MERLIN, for a detection rate of 20% consistent with previous observations 
of OH masers in other SNR’s (H03; Paper I). In addition, the D18 maser, which was not 
observed with MERLIN, is detected using the VLBA. The current VLBA observations of 
the EllA maser are consistent with the VLBA observations of C99. The current VLBA 
observations of the F24 maser (Fig. 3) indicate an extended (approximately 80 mas) angular 
structure. The F24 maser was not detected in the VLBA observations by C99, as discussed 
in §4.1. 


•^The Astronomical Image Processing System (AIPS) is documented at http://www.nrao.edu/aips/. 
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The deconvolved sizes of the MERLIN maser images are in the range 50 to 350 mas 
(150 to 1000 AU). Compact cores with deconvolved sizes of 20 mas (60 AU) are apparent 
in the VLBA images. These Endings are consistent with previous observations of OH SNR 
masers (C99; Paper I; H03). 


3.2. Line Profiles 

Tables 4 and 5 list the fitted positions and deconvolved spectral-profile properties of 
the maser lines from the MERLIN and VLBA data. All line profiles were fitted using the 
GIPSY^ software package. Figure 4 shows the full-Stokes spectra for the EllA maser from 
the MERLIN data. 

The line widths determined from the VLBA data are narrower than the line widths 
determined from the MERLIN data, a trend also observed in the IC 443 and W28 OH masers 
(H03; Paper I). (The trend of observed line width decreasing with beam size is observed for 
many maser species, e.g. OH [4765 MHz]: Palmer, Goss, & Devine 2003; H 2 CO [4830 MHz]: 
Hoffman et al. 2003b.) We assume that for OH (1720 MHz) SNR masers, since the VLBA 
data appear free from blending (see Paper I), the values for line width derived from VLBA 
data represent the intrinsic line widths of the masers. The spectra observed at the F24I and 
F24H image peaks indicate two velocity components, as summarized in Table 4. 

Figure 5 shows the residuals from the best-fit Gaussian to the Stokes / profile of the 
EllA maser from the MERLIN data. The line shapes of the W44 OH masers deviate from a 
Gaussian shape as discussed in §4.5. Also listed in the tables are the brightness temperatures 
Tb of the masers, determined using the relation Tg = (O.565A^S')/(/c0a0b), where S is the 
flux density of the maser (Wm“^Hz“^), A is the wavelength of the radiation (m), k is 
Boltzmann’s constant, and 0a and 0^ are the major and minor axes of the source (radians) 
{e.g. Rybicki & Lightman 1979). 


3.3. Zeeman Circular Polarization 

There is sufficient signal-to-noise in the current MERLIN data to measure the magnetic 
field in the maser gas using two different methods of Zeeman analysis. The two methods 
described in this section will be compared in §4.3. 


"^The Groningen Image Processing System (GIPSY) is described at http: //www. astro.rug.nl/^gipsy/. 



The Zeeman splitting in the W44 OH masers has been directly resolved as a velocity (fre¬ 
quency) difference between the right-circular polarization (RCP, using RR cross-correlations) 
and left-circular polarization (LCP, using LL cross-correlations) maser spectra. The fre¬ 
quency difference is converted to magnetic field using the relationship urcp — ^^lcp = 
where Rrl denotes the magnetic field strength measured using the frequency separation 
between the fitted RCP and LCP line profiles and Z is the Zeeman splitting coefficient 
(see Paper I for a complete discussion of the Zeeman fitting techniques). As in Paper I, 
the value Z = 0.114 kms“^mG~^ (0.654 HzpG“^; Davies 1974) is used. The observed ve¬ 
locity separation of the circularly-polarized maser components may be parametrized with 
xb = Az/^/Az/d, where Az/^ is the Zeeman frequency separation of the split substates and 
Az/d is the Doppler-broadened width of the maser line. For the W44 maser components 
Xb ~ 0.1 (Table 6). 

The following equation relates the magnetic field to the Stokes / and Stokes V line 
profiles {e.g. Roberts et al. 1993; Sarma et al. 2001; BFGT), 

dl 

V = b— + al. (1) 

The Stokes I and V profiles of the maser lines from both the MERLIN and VLBA obser¬ 
vations were fitted to equation 1. All fitting was performed with the MIRIAD^ software 
package. The al term in equation 1 represents the symmetrical component of the Stokes 
V profile. The fitted magnitude of a (approximately 0.02) is consistent with instrumental 
“leakage” of Stokes I into Stokes V] in the current data a does not represent an intrinsic 
property of the maser radiation. The values of the fitted parameter b (approximately 0.03) 
are related to the magnetic field Bq using b = ZB 0/2 (see also the discussion Paper I and 
§4.3 concerning conversions of 6 to R which involve 6). 

Table 6 lists the values of urcp — 'Wlcp; -Brl, xb, and Bq for the MERLIN observations. 
Although the W44 masers are detected in Stokes V using the VLBA, there is insufficient 
signal to noise (< 3 (t) to fit a magnetic field to the Stokes V profiles. The sign and magnitude 
of the magnetic fields fitted to the MERLIN data are consistent with VLA observations (G97). 
Figure 6 shows the good agreement between the Rrl and Bg measurements of the magnetic 
field. 


^The Multichannel Image Reconstruction, Image Analysis and Display (MIRIAD) is documented at 
http://bima.astro.umd.edu/miriad/. 
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3.4. Linear Polarization 

Table 7 lists the properties of the linear polarization of the masers observed using MER¬ 
LIN. The linear-polarization measurements have not been corrected for Faraday effects (see 
§4.4.1). Linear polarization is not detected using the VLBA with a 3-a upper limit of 10%, 
consistent with the MERLIN observations which show that the typical linear polarization 
fraction is approximately 5% to 10%. The linear polarization angles y listed in Table 7 are 
determined using the usual relation tan 2% = U/Q, where U and Q are the Stokes parameters 
formed from the cross correlations of the right- and left-circularly polarized emission. Figure 
4 shows the full-polarization prohles for the EllA maser observed with MERLIN. The dotted 
lines in Figure 4 show the fitted prohles using (1) the Stokes I line prohle summarized in 
Table 4, (2) the Stokes V prohle htted using equation 1, and (3) the htted Stokes Q and U 
prohles (Table 7). The percentage of linearly polarized intensity, q = /1, is also 

listed in Table 7. The values of y and q are in good agreement with those observed by C97 
using the VLA. 


4. Discussion 

4.1. Comparison with Previous Observations of W44 

The E and F regions of OH (1720 MHz) SNR maser emission in W44 were previ¬ 
ously observed by C99 using both MERLIN-|-Lovell and the VLBA-I-Yl. The beam size 
(290 X 165 mas at a position angle of 25°) and image noise (30 mJy beam“^) of the C99 
MERLIN images are comparable to (through slightly larger than) those of the current MER¬ 
LIN observations (240 x 141 mas at a position angle 24°, 20 mJy beam“^). The VLBA 
observations discussed by C99, however, were much shorter in duration than the current 
observations (approximately 3.5 hr on-source compared with approximately 24 hr) and the 
image noise was correspondingly higher (100 mJy beam“^compared with 25 mJy beam“^). 
The C99 VLBA images had a larger beam of 40 x 30 mas at a position angle of —10° com¬ 
pared with the beam of the current VLBA observations: 27 x 13 mas at a position angle of 
-23°. 

Using MERLIN, C99 imaged the Ell and F24 masers in W44. The angular structure 
observed in the Ell maser is similar to that shown in the current data (Fig. 2). The 
deconvolved angular size and flux density of the EllA maser are also in agreement between 
the two observations. The peak flux densities measured for the F24 maser are in agreement 
to approximately 30% (C99: 1.7 Jy beam“^; current data: 2.3 Jy beam“^) between the 
two observations. Two additional masers, E20 and F23, have been detected in the current 
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MERLIN data. Although these two masers are each about 60% weaker than the Ell and F24 
masers, the E20 and F23 masers, at their current flux densities, should have been detectable 
in the C99 data. 

Using the VLBA, C99 were only able to detect the EllA maser. In the C99 images 
the EllA maser is unresolved with a flux density of 1.2 Jy, consistent with the current 
observations of a slightly resolved source having flux density 1.1 Jy (Table 2). Based on the 
new VLBA observations, there is also a detection of the F24 maser (Fig. 3). The F24 maser 
is not detected in the VLBA data by C99, presumably due to insufhcient u-v coverage in the 
C99 observations: the 2 Jy of total flux density in the data of the F24 maser is distributed 
over three components which subtend greater than 60 mas. The C99 experiment used an 
observing schedule which divided eight hours of VLBA observing time between the W44 and 
W28 masers. The current VLBA experiment dedicated seven-hour tracks to observations of 
the W44 masers, resulting in denser u-v coverage than the C99 schedule. No polarization is 
discussed by C99 for the observations of the W44 masers. 


4.2. Maser Amplification 

The amplification process of an unsaturated maser may be parametrized by the gain 
r such that Te ^ Thge~'^, where Tb is the brightness temperature of the maser and T^g is 
the brightness temperature of the background radiation which provides the “seed” radiation 
for the stimulated emission process. In this section we discuss constraints on the agreement 
between the observed brightness temperature of the masers and the maser gain suggested 
using the collisional pump model. 

Elitzur (1976) suggests that the maximum gain of a collisionally pumped 1720-MHz 
OH maser is r ~ —20. Pavlakis & Kylafls (1996) suggest that the maximum brightness 
temperature of a collisionally pumped 1720-MHz OH maser Tg ~ 10® K, corresponding to 
a gain r ~ —20. These suggestions for the value of the gain r > —20 may be compared 
with measurements of Tb and T^g using Tb ~ Thge~'^. T^g is known because of the geometry 
of the SNR/molecular cloud interaction in W44: the molecular cloud is directly between 
the observer and the SNR [e.g. Kodaira et al. 1977) such that the masers likely amplify the 
1720-MHz synchrotron continuum radiation from the SNR. 

The brightness temperature of the background radiation at 1720 MHz may be extrap¬ 
olated from observations at 1442 MHz by Giacani et al. (1997, with resolution of 15"). For 
synchrotron radiation Tb oc {e.g. Rybicki & Lightman 1979), where a is the spectral 
index of the synchrotron radiation, measured by Giacani et al. (1997) to be a ~ —0.4. The 
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brightness temperature of the 1442-MHz synchrotron radiation in the E and F maser region 
is about 90 K, corresponding to a brightness temperature for the 1720 MHz background ra¬ 
diation of Ti,g ^ 60 K. Therefore, the gain corresponding to a maser brightness temperature 
Tb ~ 10^ K (Table 5) is r ~ —17. If the maser is saturated, this analysis yields the lower 
limit on the strength of the gain r < —17. 

For the case in which the masers amplify their own spontaneous emission, rather than 
an external background source, T^g is determined from the excitation temperature of the 
spectral line (see, e.g., Elitzur 1992). If \Tx\ > \Ti,g\ the seed provided internal to the masers 
is greater than the seed provided by the external background, in which case the gain will 
not be as strong as calculated above, indicating the upper limit r > —17. In the absence of 
hrm constraints on T^, we assume the approximate value r ~ —17. Since the estimated gain 
does not exceed the suggested maximum gain, we hud these observations to be consistent 
with suggestions based on the collisional pump model (see also §4.5.1). 


4.3. Zeeman Detections 

The magnetic held in the masers has been measured using both resolved Zeeman split¬ 
ting between the RCP and LCP line prohles (Hrl) and htting of Stokes V to the derivative 
of Stokes / {Bq). The values of Bq are consistent with the values of Hrl (Fig- 6). However, 
despite the good agreement between the derived parameters Be and Hrl, it is important 
to note that the theoretical models used to interpret the observations may not be appropri¬ 
ate. Since theoretical investigations of maser polarization are usually based on xr > 1 or 
1 (cf. Elitzur 1996), many of the theoretical models (concerning both Zeeman analysis 
and linear polarization analysis) may not be applicable to these OH SNR masers for which 
xb ~ 0.1 (§3.3; Table 6; see also Paper I). 

The sign and magnitude of the magnetic helds indicated by the current circular-polarization 
observations are consistent with previous observations of OH (1720 MHz) SNR masers. The 
sign of the htted magnetic held indicates the direction of the magnetic held projected onto 
the line of sight (although the magnitude of the Zeeman splitting may not be linearly related 
to the line-of-sight component of the held, see Paper I). All of the W44 masers are htted 
with a negative magnetic held indicating a magnetic held oriented toward the observer. The 
magnitude of the magnetic held measured using the circular-polarization data, although 
uncertain by a factor of approximately three due to uncertainties in 9 and the degree of 
saturation {e.g. Watson & Wyld 2001, hereafter WWOl), may be compared with theoretical 
suggestions about the strength of the magnetic held in shock-compressed molecular clouds. 

As discussed by BFGT, Bertoldi & McKee (1992) suggest \B\ = 0.4Ann°'®, where n is the 
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density in cm”^ and An is the shocked line width in km s“^. Frail & Mitchell (1998), using 
observations of CO (J = 3 —2), estimate the density in the post-shock to be n ~ 10^ cm“^ 
and the shocked linewidth to be An = 40 km s“^. For these values, the Bertoldi & McKee 
(1992) model yields |il| ~ 5 mG which is of the same order as the magnetic held magnitude 
estimated from the observed Zeeman splitting (^ 1 mG). 


4.4. Linear Polarization 

In this section we compare the magnetic held direction indicated by the shock front and 
the magnetic held direction indicated by the linearly-polarized maser emission. Using this 
comparison, we discuss constraints on the magnetic held geometry within the masers. 


Other Observations of Magnetic Field Orientation 

The magnetic held orientation in the SNR may be probed directly using observations 
of the synchrotron continuum radiation. Milne & Dickel (1975) have observed the linear 
polarization of the synchrotron emission from W44 at both 5 GHz and 2.7 GHz at an 
angular resolution of 4f4 (see also Milne [1987] which contains a review of the observations 
by Kundu & Velusamy [1972] and Dickel & Milne [1976]). Using a dual-frequency analysis, 
they hnd the Faraday rotation measure (RM) in the direction of the masers to be between 0 
and —10 radm“^ and the depolarization to be negligibly small. Based on these synchrotron 
observations, we conclude that the Faraday rotation and depolarization ahecting our maser 
measurements is small and need not be considered. Milne and Dickel also hnd the position 
angle of the plane-of-the-sky component of the magnetic held at the location of the masers 
to be approximately —35°. The magnetic held orientation is uniform in direction for several 
arcminutes on every side of the maser region. Based on the relative uniformity of the 
results of the arcminute-resolution synchrotron observations, we assume that the Faraday 
rotation, depolarization, and magnetic held orientation do not change appreciably over the 
< 1' angular separation of the masers. 

The coincidence of OH SNR masers with shocked molecular gas is well established {e.g. 
Frail & Mitchell 1998). As the shock forms and moves outward through the interstellar 
medium, the ambient galactic magnetic held is swept up and compressed. For example, 
Balsara et ah (2001) suggest that the direction of the magnetic held in the shocked gas lies 
in a plane parallel to the shock. There is a good agreement between the orientation of the 
magnetic held inferred using images of the shocked GO gas (Frail & Mitchell 1998) and the 
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position angle of the magnetic field measured from observations of the synchrotron radiation 
(—35°, Milne & Dickel 1975). 


4-4-2. Orientation of the Linear Polarization of the Masers 

The position angle y of the linearly-polarized emission from the OH masers in W44 
(Table 7) is parallel to the position angle of the magnetic held apparent from both the CO 
and synchrotron observations. The agreement between the orientation of the CO shock 
front and the maser position angle is shown in Figure 7. From these comparisons, it is 
apparent that the position angle of the linear polarization of the maser is parallel to the 
plane-of-the-sky component of the magnetic held in the maser gas. 

Based on theoretical investigations, Elitzur (1998) and Watson and Wyld (WWOl), 
for example, suggest that the polarization fraction q and the polarization angle y of maser 
radiation depend strongly on 6, the angle between the magnetic held and the line of sight. For 
masers viewed at an angle of 6^ < 55°, these authors suggest that the linear polarization angle 
is parallel to the plane-of-the-sky component of the magnetic held permeating the maser gas. 
For masers viewed at an angle of 6 > 55°, the linear polarization angle is suggested to be 
perpendicular to the plane-of-the-sky component of the magnetic held. The OH masers in 
W44 appear to have a linear polarization parallel to the plane-of-the-sky component of the 
magnetic held in the SNR shock, implying 6 < 55°. In contrast, the OH masers in W28 
appear to have a linear polarization perpendicular to the magnetic held of the SNR shock 
(Paper I). 

The masers are observed to have q ~ 10% (Table 7; see also Paper I; BFGT). For the 
case in which q ze 10% and in which y is parallel to the plane-of-the-sky component of the 
magnetic held, Watson and Wyld (WWOl), for example, suggest that 9 ~ 55°. Indeed, all 
observations of the linear polarization of OH (1720 MHz) SNR masers interpreted using the 
theories of WWOl and Elitzur (1998) yield 9 ze 55° (see BFGT; Paper I). 

Thus, the W44 masers and W28 masers diher in the relative orientation of the maser 
polarization and the magnetic held, for which theoretical polarization models suggest 9 > 55° 
in W28 and 9 < 55° in W44. However, since the fraction of linear polarized maser emission 
q is comparable in both the W44 and W28 OH SNR masers, 9 must be approximately 55° 
in both remnants. We consider the preponderance of the derived value 9 ^ 55° to represent 
an uncomfortable level of coincidence in the theoretical interpretations. We plan to discuss 
this coincidence further in a future paper. Also, observations with more favorable spectral 
resolution will more accurately determine line shape parameters, such as the deviation from 
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Gaussian (Watson et al. 2002) and the relative width of Stokes I and V {e.g. Vlemmings et 
ah 2002), which are suggested to depend upon the degree of saturation of the masers and 9, 
the angle between the line of sight and the magnetic held. 


4.5. Comparison with Discussion in Paper I 

4.5.1. Maser Gain Saturation 

Three separate tests indicate partial saturation of the masers. First, theoretical in¬ 
vestigations indicate that a maser must be at least partially saturated in order to possess 
linear polarization, assuming no anisotropies in the pumping [e.g. Elitzur 1998; WWOl). 
The degree of linear polarization observed in the W28 and W44 masers {q ~ 10%) is consis¬ 
tent with partial saturation (Paper I). Second, the estimated gain of the masers, r ~ —17 
(§4.2), is comparable to the maximum suggested gain for OH SNR masers, consistent with 
the conclusion that the masers are partially saturated. Third, Watson et al. (2002) suggest 
the line-shape parameters S (dehned therein) and kurtosis {e.g. Abramowitz & Stegun 1972; 
Press et al. 2001) are dependent upon the saturation state of the maser amplihcation process. 
The deviations from Gaussian line shape observed in the W44 maser prohles {e.g. Fig. 5) 
are similar to those observed for the W28 masers (5 ~ 10“^ and a positive kurtosis; note 
that Paper I inaccurately described the kurtosis as negative). This result is consistent with 
partial saturation (though it is also consistent with unsaturated masing). We intend to make 
a more quantitative analysis using the Watson et al. (2002) model in a future paper. 


4.5.2. Zeeman Analysis 

In Paper I we presented the hrst observation of Zeeman splitting in OH (1720 MHz) SNR 
masers to be resolved directly between the RR and LL line prohles. Both the MERLIN and 
VLBA observations of the W28 masers yield a resolved velocity difference between Zeeman 
components. Due to the relatively low hux density of the W44 masers, only the MERLIN 
data in this paper indicate a difference between the RR and LL line centers (Table 6). 
Nonetheless, for both the W44 and W28 observations, the Zeeman effect has been directly 
resolved between the RR and LL polarizations. 

A comparison of the magnetic held strength calculated from the RR and LL line sepa¬ 
ration (T^rl) and that calculated from the Stokes-htting method {Be) indicates, for both the 
W44 and W28 masers, that B^i^/Be ^ 1 (Eig. 6). Eor both the W28 and W44 masers, the 
ratio of line splitting to line width is xb ~ 0.1. Since Be is expected to depend upon 9 but 
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i?RL is not, we conclude that their agreement indicates that either (1) the magnitude of the 
magnetic held htted to the polarization observations does not depend as strongly on 6 as 
is currently suggested in the theoretical literature {e.g. WWOl) or (2) the observed Zeeman 
splitting in the masers {xb ~ 0.1) renders models based upon -C 1 inapplicable to these 
results. We plan to collect all polarization observations of OH SNR masers into a future 
paper which addresses the appropriateness of the theoretical interpretations. 


5. Conclusions 

Using MERLIN, we have detected all but one of the OH (1720 MHz) SNR masers in 
W28 with intensity greater than 1 Jy beam“^ as observed with the VLA by Claussen et al. 
(1997). Of the 10 MERLIN detections, we have detected two masers (approximately 20%) 
with the VLBA+Yl, consistent with other VLBI detection rates of OH (1720 MHz) SNR 
masers (H03; Paper I). Based on MERLIN data having a resolution of 200 mas and VLBA 
data with a resolution of 15 mas, the masers are measured to have deconvolved angular sizes 
of 50 to 350 mas (150 to 1000 AU) with compact cores 20 mas (60 AU) in size, consistent 
with theoretical expectation and previous observations. All of the imaging results for the OH 
(1720 MHz) SNR masers in W44 are consistent with theoretical expectation and previous 
observations. In general, OH (1720 MHz) SNR masers which are not affected appreciably by 
interstellar scattering (see Yusef-Zadeh et al. 1999a; Claussen et al. 2002) have similar sizes 
and angular morphologies (Paper I; H03). 

Zeeman splitting is observed between the right- and left-circularly polarized maser lines, 
offering a direct measurement of the magnetic field in the maser gas. The strengths of the 
magnetic helds estimated from the circular-polarization observations (~ 1 mG) are consistent 
with the magnetic held strengths suggested using models of shock-compressed molecular 
clouds. Also, the position angle of the magnetic held measured using the linear polarization 
of the masers is in good agreement with both the orientation of the shock observed using 
molecular emission [e.g. Frail & Mitchell 1998), with the magnetic held orientation expected 
from the SNR expansion [e.g. Balsara et al. 2001), and with the magnetic held orientation 
calculated from observations of the synchrotron emission [e.g. Dickel & Milne 1975). Thus, 
the linear polarization angle of the W44 OH masers appears to be parallel to the SNR 
magnetic held in contrast to the position angle of the linear polarization of the W28 OH 
masers (Paper I) which was found to be perpendicular to the magnetic held in the SNR. In 
addition, interpretation using existing theoretical models (appropriate for xb 1 although 
the masers have xb ~ 0.1) indicate that 6 ~ 55°, as has been found for all other OH-SNR- 
maser polarization observations. 
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Fig. 1.— A greyscale image of the 1442-MHz synchrotron continuum radiation from W44 
(Giacani et ah 1997). The scale at the top of the hgure is in units of mJy beam“^. The image 
has a resolution of 16" x 15". The location of the various regions of OH (1720 MHz) emission 
concentrations are indicated (after Claussen et ah 1997). The sizes of the ellipses used to 
to indicate the maser regions have no special signihcance; the regions contain between one 
and ten masers as observed with the VLA (C97). Using MERLIN we observed the masers 
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right ascension (J2000) 

Fig. 2.— Contour and greyscale image of the Ell maser region at ulsr = 44.2 km s“^ from 
the MERLIN observations. The contour levels are -4, 4, 8, 16, 32, 64, 100, and 120 times 
the image RMS noise level of 20 mJy beam“^ (no negative contours appear). The beam, 
plotted in the lower left-hand corner, is 240 x 141 mas at a position angle of 24°. 
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Fig. 3.— Contour and greyscale image of the F24 maser region at ulsr = 46.8 km s“^ from 
the VLBA observations. The contour levels are -7, 7, 14, 21, and 28 times the image RMS 
noise level of 25 mJy beam“^ (no negative contours appear). The beam, plotted in the lower 
right-hand corner, is 27 x 13 mas at a position angle of —23°. 
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Fig. 6.— A plot of the magnetic field stre 
observations (Table 6) using y/((9//(9z/)-fitt 
described in §3.3. The dashed line has a slof 
the two measurement methods. 
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Fig. 7.— Contour and greyscale image of the CO (J = 3 —2, 346 GHz) line emission at 
ulsr = 45 km s“^ from the W44 E and F maser regions observed using the James Clerk 
Maxwell Telescope (C. L. Brogan 2004, private communication). The contours are at 3, 
5, 7, and 9 times the image RMS noise level of 9 K. The beam of the telescope is 15". 
The hlled circles mark the positions of the OH masers for which the position angle of the 
linear polarization has been determined using MERLIN (Table 7); the straight lines through 
the maser markers indicate the linear polarization angle. The good agreement between the 
orientations of the shocked CO gas and the maser polarization is discussed in §4.4.1. 
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Table 1. VLB A Observation Positions 


R.A. (J2000) 

(h m s) 

Deck (J2000) 

(O , 



array tracking positions 

18 56 01.4902 

01 10 42.177 

32.2215 

24 59.361 

correlation centers 


18 56 01.2889 

01 12 46.784 

03.8015 

08 46.949 

28.1418 

29 12.076 

29.4822 

20 26.555 

36.6571 

26 33.176 









Table 2. MERLIN Maser Sizes and Flux Densities 


Feature 

^(Jy) 

angular size^ (mas) 

P.A. (deg) 

B7 

0.9 

0.24(2)x< 0.07 

173(4) 

C8A 

0.4 

0.23(5)x< 0.17 

148(21) 

C8B 

0.3 

0.34(7)x< 0.14 

171(9) 

EllA 

2.7 

0.20(l)x< 0.05 

150(3) 

EllB 

0.2 

0.37(6)x 0.15(5) 

166(12) 

Flic 

0.4 

<0.19 x<0.15 

144(18) 

E20A 

1.5 

0.35(l)x< 0.07 

155(2) 

E20B 

0.2 

0.34(6)x< 0.1 

162(8) 

F23 

0.9 

0.23(3)x< 0.11 

158(14) 

F24 

2.9 

0.19(l)x< 0.03 

165(3) 


®'the deconvolved angular size of the masers listed 
(major axis) x (minor axis) 
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Table 3. VLBA Maser Sizes and Flux Densities 


Feature 

^(Jy) 

angular size (mas) 

P.A. (deg) 

EllA 

1.2 

24(6) 

133(18) 

D18 

0.5 

< 13 

158(20) 

F24 

2.0 

19(4) 

147(9) 
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Table 4. MERLIN Fitted Maser Properties 


Feature 

R.A. (J2000) 
(h m s) 

Decl. (J2000) 

(O , 

/ 

(mjy beam“^) 

Av 

(km s“^) 

'i'LSR 
(km s“^) 

Tb 

(10® K) 

B7 

18 56 01.291 

01 12 44.80 

760 

0.65(1) 

44.63(1) 

11. 

C8A 

03.784 

08 46.97 

340 

0.93(2) 

46.16(1) 

4.7 

C8B 

03.788 

46.51 

230 

0.87(3) 

46.23(1) 

3.2 

EllA 

26.734 

29 18.15 

2410 

1.03(1) 

44.11(1) 

34. 

EllB 

26.741 

16.59 

200 

0.84(3) 

44.23(1) 

2.8 

Flic 

26.779 

17.53 

360 

0.71(1) 

44.15(1) 

5.0 

E20A 

29.684 

08.16 

940 

0.99(1) 

44.29(1) 

13. 

E20B 

29.692 

07.42 

170 

1.70(3) 

44.16(2) 

2.4 

F23 

36.633 

26 29.82 

840 

0.88(2) 

46.72(1) 

12. 

F24 

36.675 

37.14 

2390 

0.71(1) 

46.96(1) 

33. 






Table 5. 


VLBA Fitted Maser Properties 


Feature 

R.A. (J2000) 

(h m s) 

Decl. (J2000) 

(O , 

/ 

(mJy beam”^) 

Av 

(km s“^) 

t'LSR 
(km s“^) 

Tb 

(10^ K) 

EllA 

18 

56 26.7347 

01 29 18.154 

680 

0.63(3) 

43.85(1) 

0.9 

D18 


29.4822 

20 26.556 

450 

0.57(5) 

43.79(2) 

0.6 

F24AIII 


36.6741 

26 37.115 

340 

0.55(5) 

46.84(2) 

0.4 

F24AI 


36.6748 

37.136 

970 

0.59(3) 

46.79(1) 

1.3 





140 

a 

46.19(5) 

0.2 

F24AII 


36.6761 

37.100 

570 

0.48(4) 

46.77(2) 

0.7 





220 

a 

46.34(3) 

0.3 


^insufficient signal to noise to fit line width 







Table 6. 


MERLIN Zeeman Results 


Feature 

■i^RCP — '^LCP'^ 
(km/s) 

Xb 

^RL 

(mC) 

(mG) 

B7 

C 

C 

C 

-0.33±0.09 

C8A 

C 

c 

c 

d 

C8B 

c 

c 

c 

d 

EllA 

-0.084(4) 

0.08 

-0.73±0.02 

-0.75±0.02 

EllB 

C 

c 

C 

d 

EllC 

-0.066(13) 

0.09 

-0.58±0.09 

-0.57±0.08 

E20A 

-0.049(7) 

0.05 

-0.43±0.07 

-0.48±0.02 

E20B 

C 

c 

C 

d 

F23 

-0.061(11) 

0.08 

-0.53±0.09 

-0.53±0.03 

F24 

-0.061(4) 

0.09 

-0.53±0.05 

-0.56T0.02 


®'the velocity difference between the centers of the Gaus- 
sians fitted to the RCP and LCP line profiles 

’’the magnetic field strength fitted to V/{dI/dv) 

'’insufficient signal to noise to resolve Zeeman splitting 
between RCP and LCP 

'’insufficient signal to noise to observe Zeeman pattern in 
Stokes V 







Table 7. 


MERLIN Linear Polarization Results 


Feature 

Q 

(mJy beam”^) 

U 

(mJy beam”^) 

X 

(deg) 

<1 

(%) 

B7 

a 

a 

a 

a 

C8A 

a 

a 

a 

a 

C8B 

a 

a 

a 

a 

EllA 

200 

-220 

-24(1) 

12(1) 

EllB 

< 10 

< 10 

- 

< 9 

EllC 

-15 

-50 

-53(6) 

14(3) 

E20A 

25 

-30 

-26(7) 

4(1) 

E20B 

< 10 

< 10 

- 

< 8 

F23 

-60 

< 10 

0(4) 

7(3) 

F24 

70 

-150 

-32(3) 

7(1) 


’^polarization calibration unavailable 







